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ABSTRACT 

Aims. Our aim is to compare the infrared properties of big, "classical" dust grains with visual extinction in the cloud L1642. In par- 
ticular, we study the differences of grain emissivity between diffuse and dense regions in the cloud. 

Methods. The far-infrared properties of dust are based on large-scale 100 yum and 200 yum maps. Extinction through the cloud has 
been derived by using the star count method at B- and /-bands, and color excess method at J, H and K s bands. Radiative transfer cal- 
culations have been used to study the effects of increasing absorption cross-section on the far-infrared emission and dust temperature. 
Results. Dust emissivity, measured by the ratio of far-infrared optical depth to visual extinction, r(far-IR)/Av, increases with decreas- 
ing dust temperature in L1642. There is about two-fold increase of emissivity over the dust temperature range of 19 K-14 K. Radiative 
transfer calculations show that in order to explain the observed decrease of dust temperature towards the centre of LI 642 an increase 
of absorption cross-section of dust at far-IR is necessary. This temperature decrease cannot be explained solely by the attenuation 
of interstellar radiation field. Increased absorption cross-section manifests itself also as an increased emissivity. We find that, due to 
temperature effects, the apparent value of optical depth r !ipp (far-IR), derived from 100 /im and 200 /im intensities, is always lower than 
the true optical depth. 
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1. Introduction 

Recent observations have indicated that the far-IR emissivity 
(measured by the ratio r(far-IR)/Ay) increases with decreas- 
ing temperature (Cambresy et al. 2001; del Burgo et al. 120031 
Stepnik et al. [20031 Kramer et al. 120031 . Probable reasons for 
increased emissivity are coagulation of dust particles into large, 
fluffy aggregates, and formation of ices on grain surfaces in 
dense regions (Bazell & D wek [19901 P reibisch et al. [19931 
Ossenkopf & Henning [1994l Dwek [19971 Cambresy et al. |200T1 
and references therein). Variations of dust optical properties are 
visible even over a scale of galactic anticenter hemisphere, any- 
where where Ay ^ 1 mag (Cambresy et al. 2005 ). del Burgo & 
Laureijs (2005 ) used observations of the Taurus molecular cloud 
TMC-2 between 60 fim and 200 fim to separate the emission into 
a cold and a warm component. They found that far-IR emissiv- 
ity of the cold component is a few times larger than that of the 
diffuse interstellar medium, and that the change in the proper- 
ties of the big dust particles takes place at intermediate densi- 
ties of n(H2)~ 10 3 crrT 3 . They favoured grain-grain coagulation 
over gas accretion as the reason for enhanced emissivity because 
the latter process produces an increase of emissivity by a factor 
of about two at most, less than that observed by del Burgo & 
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Laureijs. Also Stepnik et al. (2003) suggest that grain-grain co- 
agulation is the reason for changes of dust properties in dense 
regions. 

For the introduction of the high galactic latitude dark cloud 
L1642 we refer to Lehtinen et al. (2004, Paper I) and references 
therein. Paper I was based on ISOPHOT observations in the far- 
infrared, and we made the following conclusions: 
Based on the 200 fim optical depth map the cloud was divided 
into three regions, called A (consists of two cloudlets, Al and 
A2), Band C (see Fig. [B- 

Re gion B, the one with the highest dust column density, Ay ~ 
8 mag, is related to a temperature minimum, * 13.8 K for a 
v 2 emissivity. 

This region is close to being gravitationally bound, while the 
others are not. 

The lower dust temperature in region B cannot be explained 
solely by the dust attenuation of the general interstellar radiation 
field, but changes in the properties of dust particles are required. 

In this paper we will compare the infrared properties of big 
grains with visual extinction. By comparing far-IR emission data 
with optical extinction data we can study the infrared and optical 
properties of dust grains. In particular, combination of reliable 
dust temperature, far-IR optical depth and visual extinction data 
enables us to study the differences of grain emissivity between 
diffuse and dense regions in the cloud. 
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Fig. 1. Visual extinction through L1642, derived using J, 77 
and K s data from 2MASS archive, smoothing the individual 
extinctions of stars with a gaussian with FWMH=4.5'. The 
IRAS sources, marked with asterisks, are from right to left: 
IRAS 04325-1419, 04327-1419, 04336-1412, 04342-1444, 
04347-1415 and 04349-1436. The FWHM sizes of the regions 
Al, A2, B and C as derived by Lehtinen et al. (Paper I) from 
200 //m optical depth data are shown as ellipses and circles. The 
contours are at 1.5, 2.5, 3.5, 4.5 and 5.5 mag 

Table 2. Error sources contributing to the pixel-to-pixel vari- 
ations of temperature, far-IR optical depth and emissivity. See 
Section 2.4 for explanations 





Error source 


Value 


i) 


< cr(I(200yum))> 


0.19 MJysr 1 


ii) 


o-(/(100yum)) 


5% 


iii) 


cr(flatfield) 


0.43 MJysr 1 


iv) 


< <r(PSF)> 


0.14 MJysr 1 


v) 


< <r(A v )> 


0.32 mag 


vi) 


< cr(zero-level(100//m))> 


0.9 MJysr 1 


vii) 


< <x(zero-level(200 fim))> 


2.6 MJysr 1 


viii) 


< cr(zero-level(A v ))> 


0.43 mag 



2. Data analysis 

For the data reduction of the ISO data we refer to Paper I. 

Extinction through L1642 has been determined via star count 
method at optical, and color excess method at near-IR wave- 
lengths. Near-IR J, 77 and K s band data have been extracted from 
the 2MASS (2 Micron All Sky Survey) archive. The two binary 
T Tauri stars near the center of the cloud have been excluded 
from the 2MASS data. 



2.1. B- and l-band extinction from starcounts 
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Fig. 2. a) Relation between visual extinction and 200 //m surface 
brightness. The data points of the regions A, B and C are marked 
with green, yellow and red, respectively. Linear fits to regions 
A/B and interclump part of the cloud are shown. We have sub- 
tracted from Ay and 7(200 yt/m) the zero-level values as fixed out- 
side the extent of L1642 (see Section 2.3.1). b) Relation between 
visual extinction and 100 pm surface brightness. Linear fits to 
the regions A/B and interclump part of the cloud are shown. The 
fit of the interclump part has been forced to go through the ori- 
gin. We have subtracted from Ay and 7(100 jt/m) the zero-level 
values as fixed outside the extent of L1642 (see Section 2.3.1) 
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Fig. 3. Relation between apparent 200 /jm optical depth and hy- 
drogen column density N(H+H.2), derived from visual extinction 
(see Sect. 3.3). Typical pixel-to-pixel 1-cr errors are shown. The 
data points of the regions A, B and C are marked with green, 
yellow and red, respectively. The values of A v and 200 fim op- 
tical depth have been derived after subtracting from Ay and in- 
tensities 7(100/um) and 7(200 jt/m) the zero-level values as fixed 
outside the extent of L1642 (see Section 2.3.1). 



Details of the starcounts are given in Appendix A. We state here 
the conclusion that extinction law in LI 642 is similar to that of 
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Table 1. Values of parameters in L1642, in other clouds and in models. For a description of the different regions of L1642 see 
Section 2.3. All the data of the cloud L1642 apply to 4.5' resolution, except columns 4 and 6 which are for 3.5' resolution. Columns 
5 and 6 give the slopes of linear fits between /(far-IR) and Ay. For a description of methods A and B see Section 3.2. Columns 2-4 
for L1642 give the values after the zero-levels of the interclump dust have been subtracted (see Section 2.3). Columns 7-10, in the 
case of method A, give the values after the background levels of each region have been subtracted (see Section 2.3). For L1642, 
the columns 2-4 and 7-10 (in the case of 'method A'), give the mean values and 1-cr standard deviations of the corresponding 
parameters within each region 











Slope of 


Slope of 




< Ay > 


< /(100 pm) > 


< /(200pm) > 


/(100 pm) vs.A v 


/(200pm) vs. Ay 


Region 


[mag] 


[MJy sr 1 ] 


[MJy sr 1 ] 


[MJy sr 1 mag -1 ] 


[MJy sr 1 mag -1 ] 


1 


2 


3 


4 


5 


6 


L1642 A 


2.1±0.6 


14±1 


34±3 


1.6+0.2 


4.2+0.4 


L1642 B 


4.3±1.1 


16±1 


52±6 


0.8+0.1 


4.5+0.3 


L1642 C 


1.5+0.2 


16±2 


35±3 


- 


- 


LI 642 interclump dust 


1.7+0.1 


12±2 


28±3 


6.7±1.6 f 


17.7+1.2 


G 300.2-16.8 ONI [1] 


1.8±0.2 




17 




9.2+2.1 


G 300.2-16.8 ON2 [1] 


1.9±0.2 




21 




11.0+2.5 


G 300.2-16.8 ON3 [1] 


2.9±0.2 




35 




12.1+2.6 
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1 dust 


< r(200yum) > 


rt200 urn) /A\, 
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no 4 i 
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Liu 11115 J 
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T \f\AT> A fmpthnrl A\ 


1 (\ 3+0 f\ 


i T+n f, 


1 3+0 ^ 


n 7+n 3 










n s+o i 


n 4+n i 




L1642B (method A) 


14.2±0.6 


4.9+2.5 


2.2±0.6 


1.2+0.3 




L1642B (method B) 


13.4±0.3 




1.9±0.3 


1.0+0.2 




T \f\AT> C impthnrl A^ 


1 a q+? f. 


U. / ITU. J 


1 1 +0 8 


n ft+n a. 




LI 642 interclump dust (method A) 


1 A T+n ^ 


a n_i_ i c\ 


Z.4+U.J 


1 3+n 3t 




\ \ f\AT) mtprrliimn Hust (mpthoH R^ 


17.1+1.0 




2.2+0.2* 


1.2+0.1* 




G 300.2-16.8 ONI [1] 


17.4 


1.9 


1.1+0.2 


0.6+0.2 




G 300.2-16.8 ON2 [1] 


17.4 


2.4 


1.3+0.2 


0.7+0.1 




G 300.2-16.8 ON3 [1] 


17.5 


4.0 


1.4+0.2 


0.7+0.1 




TPN [2] 






5.3 


2.5 




G 301.2-16.5 [3] 








0.42-0.95 




L 183 [4] 






3.4-3.8 






L 1780 [5] 


15.2 




2.6 


1.4 




Diffuse ISM [6] 






3.2 


1.7 




Diffuse ISM [7] 






2.5 


1.4 




Model [8] 








2.7 




Model [9] 






2.9 


1.5 




Model [10] 






4.7 


2.5 





* The fit has been forced to go through origin 

* Limited to visual extinction Ay > 1.2 mag 

[1] R awling s et al. ( |2004l l; [2] Lehtinen et al. ( TI998) ; [3] La ureijs et al. ( Tl996l>; [4] Ju vela et al. ( |2002l ; [5] Ridderstad et al. (2006); [6] Cambresy 
et al. < f200Tb : [7] Dwek et al. fl997l ; [8] Hildebrand ( 1983 ); [9] Li & Draine d2001bb for diffuse ISM for i? v =3.1, interpolated from their Table 6. 
Value of emissivity r(far-IR)/A v has been derived from their value of tr H (far-IR) by multiplying it with the ratio 

N(H + Hi)IA\ = 1.87 10 21 crrr 2 mag -1 (see section 3.3); [10] Desert et al. (19901 for diffuse medium. Value of <x w (far-IR) has been read from 
their Fig. 3 with about 20% accuracy. Value of emissivity r(far-IR)/Av has been derived from value of a H (far-IR) by multiplying it with the ratio 
N(H + H 2 )/A v = 1.87 10 21 cnT 2 mag" 1 (see section 3.3) 



the diffuse interstellar matter, corresponding toR y = Av/Eb-v = 
3.1 (see Fig. [2]). 



2.2. Extinction from near-IR color excess 

In order to derive extinctions from near-IR color excesses of 
stars visible through LI 642 we have applied the optimized multi- 
band technique of Lombardi & Alves ( 120011) called NICER, 
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Fig. 4. The ratio T ap p(200^m)/Av as a function of dust tem- 
perature for the regions A, B and C (green, yellow and red, re- 
spectively), as derived with the method A in Section 3.2. The 
error bars show typical pixel-to-pixel errors within each region. 
The vectors show typical direction and magnitude (1-cr standard 
deviation) of the movement of a data point within each region 
due to the inverse relation between Td us t and r(far-IR), based on 
Monte Carlo calculations (see Sect. 4.2.) 
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Fig. 6. Model clouds: Apparent dust temperature and emissivity 
(T flpp (200 fim)/Ay) as a function of multiplication factor k for the 
absorption cross-section of dust grains at far-IR (see Sect. 4.2 
for details). The dots are values of apparent temperature and 
emissivity as derived from the lOOyum and 200 //m intensities 
given by the model. The observed mean values of temperature 
and emissivity in the center of region B are marked as horizontal 
lines, shaded areas show the 1-cr uncertainties. The triangles are 
true values of emissivity 
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Fig. 5. The ratio T app (200^m)/Av as a function of dust temper- 
ature for the interclump part of the cloud, as derived with the 
method A in Section 3.2. To suppress noise, we have included 
only those datapoints for which Ay > 1.2 mag. The error bars 
show typical pixel-to-pixel errors. The straight line has a slope 
of -1.2 10~ 4 . The vector shows typical direction and magnitude 
(1-cr standard deviation) of the movement of a data point due to 
the inverse relation between Td ust and r(far-IR), based on Monte 
Carlo calculations (see Sect. 4.2.) 

which is a generalisation of the traditional color excess method 
using data of two bands only. Our data for LI 642 comprise J, H 
and K s band magnitudes from the 2MASS archive. The surface 
density of stars detected at all three bands is constant through- 
out the map, with an average of about 0.9 stars/D'. The details 
of our extinction determination are given in Appendix B. The 
resulting visual extinction map, weighted with a FWHM=4.5' 
Gaussian, is shown in Fig. Q] Region B has a maximum visual 
extinction of ~ 6.7 mag. At FWHM=3' resolution the maximum 
extinction rises to ~ 8.2 mag. Region A is seen as a local max- 



imum, but the angular resolution is not high enough to resolve 
the clumps Al and A2. In general, the correspondence between 
visual extinction on the one side and the 200 /im surface bright- 
ness and optical depth maps (see Paper I) on the other side is very 
good. There are, however, some differences. The 200 /vm surface 
brightness and optical depth of regions A and C are about equal. 
However, the maximum extinction of region C (~ 2.0 mag) is 
clearly lower than that of region A (~ 3.0 mag). Furthermore, 
the extinction over region C does not have a similar regular dis- 
tribution with a single maximum as the far-IR surface brightness 
and optical depth maps do. 

There are differences between the USNO based and 2MASS 
based extinctions. Firstly, the maximum visual extinction de- 
rived from star counts is about half of the maximum derived 
from near-IR color excesses. This can be at least partly explained 
by the lower resolution and lower penetration power of the star 
count map. Secondly, the relative values of extinction in regions 
A and C are different; in the USNO based extinction map region 
C has a larger extinction than region A, but in the 2MASS based 
map the situation is reversed. The reason for this is unknown. For 
all subsequent analyses we will use the 2MASS based extinction 
map. 

In the following data analysis of this paper, we have used 
only those datapoints for which the signal-to-noise ratio of Ay 
is greater than two. In practice, this limit sets the lowest Ay to 
about 0.5 mag. 

2.3. Zero-level of extinction and far-IR intensity 

When deriving the emissivity of dust, it is important to avoid 
biases caused by improper determination of zero level of extinc- 
tion and far-IR emission. We have determined the zero levels 
differently for the interclump dust and the dense regions. 
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2.3.1. Interclump dust 

We define the 'interclump dust' as the areas in the 200 /mi map 
which are outside the regions A/B/C, and which have Ay < 
2.0 mag. We expect that interclump dust is low density mate- 
rial located around the denser regions of LI 642, and extending 
beyond our 200/mi map. For interclump dust the zero-level has 
to be determined at locations sufficiently far away from L1642, 
outside our 200/itn map. Fortunately, in the context of another 
project, we have observed with ISO at 200/itn several positions 
in the vicinity of L1642. Here we have utilized the six posi- 
tions with the lowest 200/im intensity values, located about 2.5° 
South- West from the center of L1642 (TDT (Target Dedicated 
Time) numbers 85500550-55). 

The 100/mi intensities at the same positions have been 
derived from IRAS/ISSA data. To obtain a sufficiently low 
noise for these zero-level positions in the Ay maps, we used 
a Gaussian with FWHM=12' as a smoothing function, giving 
about 0.12 mag standard deviation per map pixel. In these six 
zero-level positions, the value of /( 1 00/im) = 3 .4 ± 0.6 MJy sr~ 1 , 
/(200/im) = 4.8 ± 0.6 MJy si-" 1 and A v = 0.14 ± 0.11 mag (the 
given uncertainty is the standard deviation of the parameter). 

2.3.2. Regions A, B and C 

We consider the regions A/B/C as separate clumps which are 
seen in excess over the general background of L1642. The zero- 
levels of regions A/B/C have been determined in the immediate 
surroundings of each of the regions. In this way, we measure 
emission and extinction which is originating solely from the re- 
gion in question. Any forground/background emission by inter- 
clump dust is in this way eliminated, at least approximately. 

For regions A and B the value of emissivity can also be de- 
rived using the slopes of /(lOO/mi) vs. /(200/mi) and /(200/im) 
vs. Ay, which are independent of any background and/or fore- 
ground zero-level subtraction. This method will be described in 
Section 3.2. 

2.4. Error estimates 

In this paper our main objective is to study the variations of the 
dust properties over the uniformly mapped L1642 cloud area. 
The different error sources contributing to the pixel-to-pixel vari- 
ations of temperature, far-IR optical depth and emissivity are 
listed in Table [2] They are: i) the intrinsic error of each pixel 
of the 200 fim intensity map, as given by the data analysis pro- 
gram PIA (ISOPHOT Interactive Analysis), ii) the pixel-to-pixel 
variation of the 100 /mi IRAS map, estimated to be 5%, iii) the 
error due to the incomplete flatfield correction of the ISO map. 
This has been estimated from the deviations of the intensities of 
adjacent pixels, in a way similar to the actual flatfield correction 
estimation (see Paper I), iv) the error due to the non-matching 
point spread functions of the 100 /mi IRAS map and the 200 /mi 
intensity map which has been convolved to the IRAS resolu- 
tion. This error has been estimated by convolving the ISO map 
to 5.0' resolution instead of the 4.5' resolution. As expected, 
the largest differences, about 1.6 MJy sr" 1 , are at the positions 
of the largest gradients at region B. This error source has been 
applied to the 200 fim map only, v) the error of the Ay map as 
given by the NICER program, vi) the error related to uncertainty 
of zero-level subtraction of 100 fim intensity (see Section 2.3), 
vii) the error related to uncertainty of zero-level subtraction of 
200 /im intensity, viii) the error related to uncertainty of zero- 
level subtraction of Ay. We have used the standard deviations of 



the zero-level values of each region as estimates of the last three 
error components. 

3. Results 

3.1. Relation between visual extinction and far-IR surface 
brightness 

The relations between 200 fim and 100 /mi surface brightness 
and Ay are shown in Fig. [2] There is a clear difference between 
the overall relations; the 1(200 /mi) vs. Ay relation shows less 
scatter, and the change of slope at Ay « 2 - 3 mag is less pro- 
nounced than for 1(100 //m) vs. Ay. 

The parameters of linear fits to the data points 1(1 00/im) vs. 
Ay and 1(200 /im) vs. Ay of different regions are given in TableQ] 
For region C, there is no clear correlation between far-IR in- 
tensity and Ay, mainly because the Ay-range in this region is 
rather small (see Fig. [TJ. It is notable that the highest values of 
1(1 00/im) are located in region C where the visual extinction, 
Ay * 1.5 mag, is lower than in regions A and B (see Fig. |2j>). 
This is due to the higher temperature in region C (see Figs. 3 and 
4 of Paper I). 

The general distribution of the data points is more linear at 
200/im than at 1 00/im. This can be explained by the fact that 
the temperature drop in region B influences the 100 /mi inten- 
sity more strongly than the 200 fim intensity because the 100 fim 
intensity is on the short-/! side of the emission maximum. 

3.2. Determination of apparent emissivity t 3PP (200 fim)/A v 

When investigating the possible variations of the dust emissivity 
in LI 642 it is important to recognize that the optical depth as de- 
rived from the far-infrared observations at 1 00/im and 200 /mi, 
i.e. 

T(200//m) = 7 v (200//m)/B v (r dust ) (1) 

is in general not the true optical depth of dust integrated over 
the line of sight. We designate the two quantities by apparent 
optical depth, T app (200/im), and true optical depth, r true (200/im) 
(see Section 4.3 for a detailed discussion). 

The relation between apparent 200 /im optical depth and Ay 
is shown in Fig. [3] At 1 00/im the relation is exactly the same 
except that the T app (100/mi) values are higher by a factor of 
(200/ 100) 2 = 4. The Pearson correlation coefficients for the re- 
lation are 0.85 and 0.91 for the regions A and B, respectively. 

We have determined the apparent emissivities using two dif- 
ferent methods. 

Method A. In this method we first subtract from the 100 /im and 
200 /im intensities and visual extinction the zero level values cor- 
responding to each region (see Section 2.3). Then the tempera- 
ture and 200/im optical depth are derived, on a pixel-by-pixel 
basis by using Eq. 1 . Within the different regions we derive the 
following mean values of apparent emissivity: 
region A: 1.3±0.5 10~ 4 mag _1 
region B: 2.2+0.6 lO^mag" 1 
region C: 1.1+0.8 lO^mag" 1 
interclump dust (Ay >1.2mag): 2.4+0.5 10 _4 mag _1 
The Ay values of interclump dust have been limited to values 
greater than 1.2 mag in order to suppress noise. The given error 
of emissivity is the l-cr standard deviation of individual data- 
points within each region. The values of emissivity at 100 /im are 
a factor of four ((200/100) 2 ) higher than the values at 200 //m. 
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Method B. From the slope of the 7(200 //m) vs. 7(100 jum) re- 
lation we derive the dust temperature 7d us t, not affected by any 
constant foreground and/or background emission components. 
This temperature represents a typical temperature in each region, 
and is given in column 7 of TableQ] Then, the 7(200 yum) vs. Ay 
relation (Fig.|2^) gives the slope between these quantities, called 
k. The apparent emissivity e app is then e app = k/B(Ti ust ). This 
method has been applied to regions A and B and interclump dust 
for which we can make a reliable linear fit of 7(200 fim) vs. Ay. 
We derive the following mean values of apparent emissivity: 
region A: 0.8+0.1 10" 4 mag" 1 
regionB: 1.9+0.3 10~ 4 mag" 1 
interclump dust: 2.2+0.2 10 _4 mag _1 

These values are about a factor of 0.6-0.9 lower than those de- 
rived with method A above. More importantly, both of these in- 
dependent methods give the same relative difference between re- 
gions A, B, and interclump dust. 

3.3. Apparent absorption cross section per H-atom cr" pp (A) 

In this section we shall derive the ratio between the far-infrared 
optical depth r app (/l) and the hydrogen column density, i.e. the 
apparent average absorption cross section per H-nucleus, 

<*> (A) - im {) 

We can use the near-IR extinction for estimation of the to- 
tal hydrogen column density, N(H + 772). We adopt the value 
N(H + H 2 )/E(B - V ) = 5 .8 10 21 cirT 2 mag" 1 valid for diffuse 
clouds (Bohlin et al. [19781 together with A V /E(B - V) = 3.1 
to obtain N(H + 77 2 )/A v = 1.87 10 21 cm" 2 mag" 1 . The rela- 
tion between N{H + 772) and r app (200 / um) is shown in Fig. [3] 
Within the different regions we derive the following mean val- 

ues of cr" p (200/im), calculated with the method A above (see 
Section 3.2): 

region A: 0.7±0.3 10~ 25 cm 2 H-atom 1 

region B: 1.2+0.5 10" 25 cm 2 H-atom -1 

region C: 0.6+0.4 10~ 25 cm 2 H-atoirT 1 

interclump dust: 1.3+0.3 10~ 25 cm 2 H-atorrr 1 

The given error of cr™ p (A) is 1-cr standard deviation of individ- 

ual datapoints within each region. The values of cr" p (100/vm) 

are four times as high as the values of cr^ p (200jum). 

3.4. T app (200 fim)/A v as a function ofT dmt 

Figure [4] shows the emissivity T app (200jum)/Ay as a function of 
dust temperature for the regions A/B/C, as derived with method 
A (see Section 3.2). The emissivity appears to increase as the 
temperature decreases. Taken from Table [2] there is two-fold 
increase of emissivity in the temperature range ~19-14K, the 
range determined by the mean temperatures of regions C and B. 
At the very lowest temperature values, 12-14K, there is indica- 
tion of a further increase of emissivity, but the number of data 
points is not sufficient to confirm that. 

The emissivity as a function of dust temperature for the inter- 
clump dust is shown in Fig. [5] The emissivity appears to increase 
by a factor of about 1.5 when the dust temperature decreases 
from 17 K to 15.5 K. 

In order to confirm the validity of the relation between emis- 
sivity and dust temperature, we must exclude the possibity that 
the relation is due to uncertainties in the derived dust tempera- 
ture: because the optical depth depends inversely on dust tem- 



perature (t oc 7/7?dust) the x- and y-axes in Figs.|4]and|5]are not 
independent. Any overestimation of temperature causes under- 
estimation of optical depth and vice versa and as a consequence 
a relation like that in Figs.|4]or|5]may occur. We have estimated 
the uncertainty of the emissivity via Monte Carlo calculations, 
considering the eight error sources as listed in Table [2] For de- 
tails see Appendix C. We conclude that the increase of emissiv- 
ity with decreasing Td us t within region B (yellow points in Fig.|4| 
can be partly caused by errors in the Td us t estimates. However, 
the increased emissivity in region B relative to regions A and 
C cannot be explained this way. It is found to be a real effect 
caused by the change of dust properties. On the other hand, in 
the case of interclump dust the range of temperature in Fig. [5] is 
rather small, and thus the change of emissivity is only slightly 
larger than the derived uncertainty. 



4. Discussion 

4.1. Dust emissivity variations 

The relatively large area of 1.6d° which we have mapped in 
and around LI 642 contains both dense and more diffuse areas 
with different dust populations. The 200 /vm and 100 jum emis- 
sion, dominated by the large particles, and the 60yum, 25 /im 
and 12/vm emissions, dominated by smaller particles, peak at 
separate positions in the cloud (see e.g. Laureijs et al. 1987). 
Compared with the situation where a sample of distinct clouds 
have been studied to cover a range of different dust properties 
(see e.g. del Burgo et al. 2003) we have the observational advan- 
tage that a more accurate (relative) calibration across our single, 
uniformly mapped area can be achieved. Thus, the variations of 
the dust properties over the mapped area are more accurately de- 
termined. 

We have found in Sect. 3.4 that the emissivity of the dust, 
T app (200 / Km)/Ay, is not constant but decreases as a function of 
dust temperature, Tt us t- Part of this variation is caused by the fact 
that the x- and y-axes are not independent. Our error analysis in 
Appendix C has shown, however, that on the large scale over the 
whole cloud area these tendencies are real. 

Based on the temperature only, the emissivity of interclump 
dust (Tdust ~ 16.3 - 17.1 K) is expected to be at about the same 
level as the emissivity of region A (Td us t ~ 15.9 - 16.3 K). 
However, it is about two times higher, at the same level as the 
emissivity in region B. This can be explained in terms of differ- 
ent contributions by large grains in different regions of LI 642. 
The maps of Laureijs et al. ( 1987 ) show that 25 /im emission has 
local maxima at regions A and C, while 60 fim emission has a lo- 
cal maximum at region C. This indicates that very small grains 
(VSGs) have a larger abundance in regions A and C as compared 
to other regions. One possible explanation for the lower far-IR 
emissivity of big grains as observed in these regions A and C 
could be the following: the big grains and VSGs are part of the 
same grain size distribution; decrease in the number of VSGs 
means an increase in the number of big grains, e.g. via coagula- 
tion or mantle formation. 

There is, in addition, the question whether the 100 /mi emis- 
sion in regions A and C contains a contribution from tran- 
siently heated very small grains, which would make the derived 
temperature an ill-defined quantity. This has been discussed in 
Appendix D, and the conclusion is that there is not evidence for 
substantial emission at 100 /mi from transiently heated VSGs in 
regions A or C. 
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4.2. Dust emissivity versus temperature: modelling 

In Paper I we showed that the decrease of dust temperature to- 
wards the centre of region B in L1642 can not be explained 
solely by the attenuation of the radiation field. In the denser 
parts of a cloud the dust grains are expected to stick together, 
forming grains with a bigger cross-section. We have used the 
radiative transfer program of Juvela & Padoan (2003 ) to test if 
an increased value of absorption cross-section (that is, increased 
emissivity) of dust grains can explain the observed decrease of 
dust temperature in LI 642. 

The parameters of the model cloud which are based on obser- 
vations are the cloud size, radial density distribution, and visual 
extinction Ay at the cloud center. Our model is tailored for re- 
gion B in LI 642. The FWHM size of r(200/mi) of the model is 
equal to the observed one. The density is constant up to 0.2 times 
the cloud radius, outside of which it has the form p(r) oc r~ L5 . 
This density distribution is in agreement with the CO observa- 
tions (Russeil et al. 2002 ) and the far-IR observations (Paper I). 
The cloud is heated by the Solar neighbourhood interstellar ra- 
diation field (Mathis et al. 119831 ). The grain properties used in 
the model are those of Li & Draine (2001b!). We have calculated 
a set of five models, multiplying the values of absorption cross- 
sections given by Li & Draine by a factor of £=1.0, 1.5, 2.0, 3.0, 
4.0, and 5.0, at all wavelengths longer than an arbitrarily chosen 
threshold wavelength of 50//m. 

The starting point of the model is the observed maximum 
visual extinction of Ay = 6.7 mag at 4.5' resolution. Firstly, in 
the case of £=1.0, the column density of dust is modified until 
the model gives a value of Ay which is equal to the observed 
one; Ay(model) is derived by integrating the extinction coeffi- 
cient along the line of sight through the model cloud. Then, in 
subsequent calculations, with increasing values of k, the true col- 
umn density of dust in the model is kept the same as in the case 
£=1.0. The value of visual extinction given by the model is the 
same for all values of k because column density of dust is con- 
stant, and the properties of dust particles are changed only at 
A > 50//m. 

The model gives the 100/mi and 200//m intensities at differ- 
ent radial distances from the cloud center. The intensities are 
convolved to 4.5' resolution. The optical depth derived from 
100/mi and 200/mi intensities given by the model cloud is an ap- 
parent optical depth, because it is obtained by mimiking the ac- 
tual ISO observations. The true optical depth, on the other side, 
is derived by integrating the absorption coefficient along the line 
of sight through the model cloud. We then derive maps of color 
temperature and apparent optical depth r app (200 /mi). The values 
of temperature and emissivity, T app (200//m)/Ay, at the centre of 
the model cloud are shown in Fig. [6] as dots for different val- 
ues of k. The observed dust temperature and emissivity at the 
center of region B are marked as horizontal lines. The shaded 
areas show l-cr uncertainty of the observed values, derived by 
assuming 20% absolute uncertainty for the 100/zm and 200/mi 
intensities. Clearly, values of absorption cross-sections at far-IR 
enhanced by a factor of about 2-3 are required to explain both 
the observed temperature and apparent emissivity towards the 
center of the cloud. The true values of emissivity, at 4.5' resolu- 
tion, are shown as triangles. 

In the case k - 2 the model predicts intensities at the 
cloud center of about 15 MJy sr~' and 64 MJy sr -1 at 100//m and 
200/mi, respectively, at 4.5' resolution. These values compare 
well with the observed ones, ~14MJy sr~' and ~60MJy sr" 1 . 



4.3. The difference between true and apparent emissivity 

As discussed in the previous section, the apparent value of emis- 
sivity is always smaller than its true value. The underlying rea- 
son is that the color temperature is biased towards too high val- 
ues, leading to an apparent optical depth which is too small. This 
bias is caused by the following facts: i) there is a temperature 
gradient towards the cloud center with lower temperatures at the 
centre than at the edge (see eg. Zucconi et al. 2001 and Galli 
et al. ( 120021 ) for a theoretical discussion), ii) graphite particles 
are warmer than silicate particles which furthermore extends the 
temperature range of dust particles, iii) the very small grains, 
which have a larger temperature fluctuation range than the big 
grains, have a noticeable contribution at 100 /mi. The net re- 
sult of the above effects is that the emission is to some degree 
dominated by the warmer dust particles; they contribute a larger 
fraction of the total emission than what should be accounted for 
by their relative mass. This is especially true at 100 /mi. The 
colour temperature derived from 100/mi and 200 /mi observa- 
tions is therefore larger than the proper, mass-weighted temper- 
ature. Thus, the apparent optical depth T app (200/mi) is always 
smaller than the true optical depth T true (200/mi). Consequently, 
the apparent emissivity is always smaller than the true emissiv- 
ity. 

Irrespective of the dust temperature, an increase of absorp- 
tion cross-section of dust particles produces an increase of opti- 
cal depth by the same amount. In Fig.[6j the apparent emissivity 
is seen to increase only by a factor of about 3.5 when k increases 
from one to five. This is due to the fact that the above mentioned 
effects are more pronounced when the temperature is lower. 

4.4. Comparison with grain models and other clouds 

The observed two-fold increase of emissivity in L1642 is sim- 
ilar to that observed in other regions; increase at 850/mi by a 
factor of ~4 over a temperature range ~20K-12K in the dark 
cloud IC 5 146 (Kramer et al. 120031 and by a factor of ~2 in B68 
(Bianchi et al. 2003); increase in the FIR by a factor of ~3.4 in 
a dense filament in Taurus (Stepnik et al. 120031 ), and increase by 
a factor of >4 in high-latitude clouds (del Burgo et al. 120031 ). 

TableQ]gives a compilation of some observed and theoretical 
properties of dust. The theoretical values for r(200/mi)/Ay and 
cr H (200 /mi) in Table Q] are "true" values in our nomenclature. 
Because the "apparent" values of cr H (far-IR) and emissivity are 
always smaller than the true values (see Sect. 4.3), a comparison 
with the model results is not straightforward. In fact, one should 
use a certain model to derive intensity values at the two far-IR 
wavelengths, 100 /mi and 200 /mi, and use them to derive the 
apparent optical depth and emissivity as we have done in Section 
4.2. 

The value of apparent cr^ (200/mi) of region B is slightly 
lower than the true value obtained for diffuse ISM by Li & 
Draine (2001a) for their carbonaceous-silicate grain model in 
the case /?v=3.1. By using Fig. [6] we can transform the value 

of cr" p (200 //m) of region B into its true value. By assuming that 
the multiplication factor for emissivity k = 2.5, the ratio of true 
emissivity to apparent emissivity for the model is about 2.5. The 

value of crjjj e (200/mi) of region B would then be about 3.0 10~ 25 . 
This is higher than any value for diffuse ISM in TableQ] 

The regions A and C have < Ay >=1. 5-2.1 mag, slightly 
higher than the extinction through typical diffuse dust. In these 

regions the value of apparent cr^J p (200 //m) (as measured with 
method A) is clearly less than the true value obtained for dif- 
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fuse ISM by Li & Draine. However, we can not compare the 
estimated values of true <x^ e (200 //m) of regions A and C to the 
true value of the model by Li & Draine because we lack radiative 
transfer calculations for regions A and C. 

Mantle accretion on dust grains can increase emissivity by a 
factor of about two at most (Preibisch et al. 1993), while grain- 
grain coagulation (Bazell & Dwek [19901 can increase the emis- 
sivity by a larger amount. Thus, either of these models can ex- 
plain the observed increase of emissivity by a factor of two be- 
tween regions B and C. 

5. Conclusions 

Based on a large scale 200 /vm ISOPHOT mapping of the cloud 
Lynds 1642, combined with 100 fim IRAS/ISSA map, and an 
extinction map we draw the following main conclusions: 

- We have found that one has to make a distinction between 
apparent values of parameters such as temperature, optical 
depth and emissivity, and their true values. The apparent val- 
ues are those which have been derived from the intensity 
maps at 100 /mi and 200 /urn. We find that the apparent op- 
tical depth and emissivity are always smaller than their true 
values because the apparent temperature is biased towards 
larger values 

- The apparent emissivity T app (far-IR)/Av increases with de- 
creasing dust temperature. There is about two-fold increase 
of apparent emissivity when dust temperature decreases 
from 19 K to 14 K 

- Radiative transfer modelling shows that an increase of ab- 
sorption cross section of dust at far-IR is capable of explain- 
ing two phenomena observed in L1642: i) decrease of dust 
temperature towards the centre of the dense region B by an 
amount that is more than what can be explained solely by 
the attenuation of interstellar radiation field, ii) increase of 
apparent emissivity towards the colder regions. We find that 
an increase of the far-IR cross sections by a factor of about 
2-3 is required to explain both the temperature drop and the 
observed emissivity increase in region B 
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Fig. .2. Relation between star count based B- and /-band extinc- 
tion. The data points of region A, B and C are plotted in green, 
yellow and red, respectively. The slopes of the dashed lines cor- 
respond to the ratios of total to selective extinction Ry = 3.1 
and R v - 5.0. The solid line is a linear fit for A B < 3.5 mag. 
The triangles are mean values of extinctions in 0.5 mag B-band 
intervals up to 4.0 mag, and the dotted line is a linear fit of type 
y — bx to these mean values for Ab < 3.5 mag. Some examples 
of 1-cr errors are shown 
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Fig. .3. The intensities measured at a position inside region C, 
at 100 /mi, 120 /mi and 200 /urn. The solid curve is a modified 
black-body with a v 2 emissivity, fitted to the 120 /mi and 200 /mi 
data points. The thick vertical errors at 100 /mi show the 1-cr 
uncertainty of the fitted black-body, assuming 10% error for the 
120 /mi and 200 /mi intensities. The thin vertical bars show the 
estimated errors of the observed intensities, 10% for the ISO data 
at 120 /mi and 200 /mi, and 20% for the IRAS data at 100 /mi. 
The horizontal bars show the bandwidth of each filter 



Appendix A: B- and /-band extinction from 
starcounts 

Data for star counting at B- and /-bands have been extracted 
from the USNO-B1.0 archive (Monet et al. 2003). We have used 
data from the second epoch SERC-J (Bj-band) and SERC-I (/- 
band) surveys. For the effective wavelengths of the B and / bands 
we have used 470 nm and 8 10 nm, respectively (see Evans 1989). 



We were forced to use a reference area relatively close to the 
L1642 (at R.A.(J2000)=4 h 35 m 00 s , Dec.(J2000)=-16°00'00", 
/=212.90°, £=-37.26°, a circular region of 1 "radius) in order to 
have the reference stars on the same plate as most of the field 
stars, and thus with the same magnitude calibration. The value 
of E(B-V) at the reference position is about 0. 12 mag according 
to the dust map of Schlegel et al. ( 1998), which corresponds to 
Ay ~ 0.4 mag. Star counting was performed in circles having a 
radius of 3'. Extinctions were derived from the \ogN{m) vs. m 
curves in the cloud and reference areas using the method of Wolf 
(I1923I I. The pixel size of the extinction map was chosen to be 3'. 

Figure [Tl shows the B-band extinction map, after smoothing 
the original map with a Gaussian with FWHM=4.5'. The max- 
imum extinction in the smoothed map, A(B) «3.0mag, occurs 
close to the center of region B. In the center of region C the 
extinction is about 1.8 mag, while in regions Al and A2 the ex- 
tinctions are about 0.8 mag and 1.0 mag, respectively. 

The relation between the B- and /-band extinctions is shown 
in Fig. [2] The extinction ratios at B- and /-bands (470 nm and 
810nm) for the cases Ry - 3.1 (diffuse dust) and Ry - 5.0 
("outer-cloud" dust) have been adopted from Table 1 of Mathis 
(1 19901 1, and are shown in Fig. [2] as dashed lines. These values 
of Ry correspond to slopes of 0.48 and 0.58 in Fig. [2] The solid 
line is a fit of type y — a + bx for all the data points which have 
Ab < 3.5 mag, weighting by the error of extinction, giving a 
slope 0.54±0.02. 

The triangles are mean values of extinctions within 0.5 mag 
B-band intervals from Ab=0.0-0.5 mag up to 3.5-4.0mag, and 
the dotted line is a linear fit of type y — bx to these mean values 
for Ab < 3.5 mag. The derived slope is 0.49+0.02. We may thus 
conclude that the extinction law in LI 642 is similar to that of 
diffuse interstellar matter. 

There is indication in Fig.[2]that region C has a lower value 
of Ry than region B. In addition, at Ab £ 3 mag the relation 
indicates a higher value for Ry inside region B. The number 
and accuracy of data points is not sufficient for firm conclusions, 
though. 

Appendix B: Extinction from near-IR color excesses 

As reference fields, representing areas with "no extinc- 
tion", we have selected two 40' diameter circular areas 
located closely at the same galactic latitude as LI 642, 
at R.A.(J2000)=4 h 32 m 6 s , Dec.(J2000)=-12°50'27", and 
R. A. ( J2000) =4 h 4 1 m 22 s , Dec. ( J2000) =- 1 6 °4 1 '44" . 

The observed colors of stars in the reference areas have the 
mean values and standard deviations of (J - Z/)=0.52+0.17 and 
(H - K) =0.19 ±0.19. They will be used as estimates of the intrin- 
sic colors ( J - H)o, (H - K)o and their variances. For the ratio of 
visual extinction to color excess we have used Ay/E(J - H) — 
8.86 an d A V /E (H - K s ) = 15.98, which correspond to R v = 3.1 
(Mathis [T990l . 

The pixel size in the extinction map was chosen to be the 
same as the pixel size of ISO C200 detector, i.e. 1.5'. The extinc- 
tion value for each map pixel was derived by using a Gaussian 
as a weighting function for the extinction values of individual 
stars, and applying the sigma-clipping smoothing technique of 
Lombardi & Alves (2001). The size of the Gaussian function 
depends on whether we compare visual extinction with 200/mi 
ISO map or with maps from IRAS data. In the case of ISO 
map we use a Gaussian with FWHM=3.5', limited by the ne- 
cessity of having a sufficient signal-to-noise ratio for the ex- 
tinction. The point-spread profile of ISO at 200//m can be de- 
scribed by a Gaussian function having FWHM=86". Thus, the 
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ISO map has been convolved with a Gaussian having FWHM= 
V(210" 2 - 86" 2 ) = 192" when comparing with the visual ex- 
tinction map having FWHM=3.5' resolution. 

The temperature and optical depth maps have been made 
by utilizing IRAS data which have a FWHM=4.5' resolution, 
measured from intensity profiles of point sources in the IRAS 
map. For these maps the 200//m ISO data have been con- 
volved with a Gaussian having FWHM= V(270" 2 - 86" 2 ) = 
256". Consequently, we use the visual extinction map with 
FWHM=4.5' resolution in the context of temperature and op- 
tical depth maps. 

There are two sources of error for the extinction map; vari- 
ance of the intrinsic colors (J - H)q and (H - K)o, and vari- 
ance of the observed magnitudes of the programme stars. In the 
Ay map made with a FWHM=4.5' Gaussian the former source 
dominates with a ~ 0.29 mag l-cr error per pixel, while the lat- 
ter source contributes typically ~ 0.08 mag l-cr error per pixel, 
giving a typical total error of ~ 0.32 mag per pixel. 

Appendix C: Error analysis of the emissivity 

We have estimated the uncertainty of the 200 /mi emissivity via 
Monte Carlo calculations, considering the eight error sources as 
listed in Table [2] For the emissivity estimation, the largest er- 
ror source is the error of Ay because it is a large fraction of 
the extinction (~ 10-75% for Ay=6-1 mag). On the other hand, 
the error of Ay does not cause any inverse emissivity vs. Td us t 
relation. Typical errors of emissivity and temperature, based on 
Monte Carlo calculations, are plotted in Fig. [4] for each of the 
region A/B/C, and for interclump dust in Fig. [5] 

We have used the Monte Carlo method to study into which 
direction and how much a datapoint moves in a Td us t vs. r(far- 
IR) plot due to pixel-to-pixel errors in the 100//m and 200 /mi 
intensity maps. Let us first consider the non-realistic case where 
only one of the intensity maps has errors. In this case there is 
always a negative correlation between Td us t and r(far-IR), just 
as expected. In the case where both of the intensity maps have 
errors, also positive correlation between Td us t and r(far-IR) is 
possible. This happens when 100 /mi and 200 /mi intensities ei- 
ther increase or decrease simultaneously. For example, if both 
of them increase by a suitable amount, the derived temperature 
may increase just slightly, and then the optical depth increases 
too because the relative increase of intensity may be greater than 
the influence of the increase of temperature. However, a positive 
correlation between Td us t and r(far-IR) happens only when the 
change of temperature is minuscule. When taking into account 
all the errors of Table|2j we find that the slope of the inverse rela- 
tion between Td us t and emissivity has a value of -0.5+0.2 10~ 4 , 
-1.1 ±0.3 10 4 and -0.4±0.3 10 4 in region A, B and C, respec- 
tively. Thus, the value of the slope increases as the temperature 
decreases. For interclump dust the predicted slope has a value of 
-1.0 ± 0.4 10~ 4 , close to the observed value of -1.2 10~ 4 . 

The vectors in Figs. [4] and [5] show the direction and mag- 
nitude (l-cr standard deviation) of the movement of datapoints 
caused by the inverse relation between Td us t and emissivity, 
when taking all the error sources of Table [2] into accout. 

Appendix D: Existence of emission from 
non-equilibrium very small grains 

Because the abundance of very small grains (VSGs) is enhanced 
in regions A and C, we have to investigate whether there is 
any significant emission from these transiently heated grains at 



100 //m. We assume that 100 /mi is the longest wavelength which 
may be affected by non-equilibrium grains. This is supported by 
the observations of L1780 (Ridderstad et al. 120061) . which is a 
high latitude cloud similar to LI 642, and the translucent/cirrus 
cloud as studied by Laureijs et al. (2000). In both these clouds, 
the intensities at 100 pun, 120 /mi, 150 /mi and 200 /mi are all 
well fitted with a single black-body function, even in the po- 
sitions where 12 /mi and 60 /mi emission have their maximum 
values. On the other hand, at 80//m there is a significant excess 
over the black-body emission from thermal equilibrium grains. 

In Fig. [3J we plot 120 /mi and 200 /mi intensities at the po- 
sition in region C, together with a modified black-body function 
fitted to them, and an estimated error of the black-body func- 
tion at 100 /mi. When taking the errors into account, there is no 
evidence for an excess emission at 100//m. We believe that any 
excess emission in region A is even smaller, because region A is 
colder than region C. 

The 120 //m and 200 //m points in Fig. [3J are from ISO 
data, the 100 /mi point is from IRAS data (scaled to DIRBE 
calibration as in Lehtinen et al. 120041 ) . The ON-source po- 
sition in region C is at coordinates R.A.(J2000)=4 h 36 m 7.2 s , 
Dec.(J2000)=-14°36'3". From the ON-source (ISO TDT 
numbers 83901895, 83901887) intensities we have subtracted 
the mean values of intensities at two reference positions, 
located at coordinates R.A.(J2000)=4 h 33 m 17.3 s , Dec.(J2000)=- 
15°58'12"(ISO TDT numbers 83901881, 83901882, 83901890, 
83901889, 83901885, 83901893) and R.A.(J2000)=4 h 28 m 59.4 s , 
Dec.(J2000)=-15°00'll"(ISO TDT numbers 83901896, 
83901888). 



